Waves attractors in rotating fluids: a paradigm for ill-posed Cauchy problems 
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In the limit of low viscosity, we show that the amplitude of the modes of oscillation of a rotating fluid, 
namely inertial modes, concentrate along an attractor formed by a periodic orbit of characteristics of 
the underlying hyperbolic Poincare equation. The dynamics of characteristics is used to elaborate 
a scenario for the asymptotic behaviour of the eigenmodes and eigenspectrum in the physically 
relevant regime of very low viscosities which are out of reach numerically. This problem offers a 
canonical ill-posed Cauchy problem which has applications in other fields. 
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Rotating fluids encompass all fluids whose motions are 
dominated by the Coriolis force. These flows play an 
important role in astrophysics or geophysics where the 
large size of the bodies makes the Coriolis force a promi- 
nent force. Some engineering problems like the stability 
of artificial satellites also require the study of rotating 
fluids because of their liquid-filled tanks M . This latter 
problem is related to the existence of waves specific to 
rotating fluids, namely inertial waves, which easily res- 
onate. These waves play also an important part in the 
oscillation properties of large bodies like the atmosphere, 
the oceans, the liquid core of the Earth ||, rapidly rotat- 
ing stars II or neutron stars H . As such, they have been 
considered since the work of Poincare on the stability of 
figures of equilibrium of rotating masses ^). Pressure 
perturbations of inertial modes for inviscid fluids obey 
the Poincare equation (PE) (christened by Cartan ||) 
which reads AP - (2tt/uj)- 2 d 2 P/dz 2 = where Qe z is 
the angular velocity of the fluid and lo is the frequency 
of the oscillation. Since ui < 2f2 M, the PE is hyper- 
bolic (energy propagates along characteristics) and since 
its solutions must meet boundary conditions, the prob- 
lem is ill-posed mathematically. Although some smooth 
solutions exist (for instance for a fluid contained in a 
full sphere or a cylinder), one should expect singular so- 
lutions in the general case. These latter solutions have 
been made explicit only recently thanks to numerical sim- 
ulations which include viscosity to regularize the singu- 
larities and let this parameter be very small as in real 
systems [§§. 

In this letter we wish to present a scenario, based on 
analytical and numerical results, for the asymptotic be- 
haviour of inertial modes at small viscosities. We use the 
case of a spherical shell as a container, which is relevant 
for astrophysical or geophysical problems, but it will be 
clear that this case is general. We will only sketch the 
main results, more details can be found in While 
the fluid mechanical problem is of much interest by it- 
self, it opens new perspectives in the theory of Partial 



Differential Equations (PDE) and also offers a toy model 
for some (very involved) problems of General Relativity 
which we shall present briefly. 

The model we use is a spherical shell whose inner radius 
is rj R and outer radius R (r/ < 1). The fluid is assumed 
incompressible with a kinematic viscosity v. We write 
the linearized equations of motion for small amplitude 
perturbations for the velocity u in a frame corotating 
with the fluid; momentum and mass conservation imply: 

du = „ . 

— + e z x u = -Vp + EAu, V-u = (1) 

when dimensionless variables are used; (20) —1 is the time 
scale and E = v/2Q,R 2 the Ekman number. When E is 
set to zero and u is eliminated, one obtains the Poincare 
equation. In nature E <C 1 and one is tempted to 
use boundary layer theory and singular perturbations to 
solve (0). However, this is feasible only when regular so- 
lutions exist for E = 0; this is the case when the container 
is a full sphere but not when the container is a spheri- 
cal shell. Indeed, numerical solutions of the eigenvalue 
problem issued from (^), where solutions of the form 
u(r)e xt are searched for (with — 1 < uj = Im{\) < 1), 
yield eigenmodes of the kind shown in Fig. 0. In this 
figure we see that the amplitude of the mode is all con- 
centrated along a periodic orbit of characteristics of the 
PE; we found this property to be quite general, after ex- 
tensive numerical exploration of least-damped modes of 
(0) i an d will now explain its origin and consequences 
on the asymptotic spectrum of inertial modes. For this 
purpose we will use axisymmetric modes since the az- 
imuthal dependence of solutions can always be separated 
out because of the axial symmetry of the problem. 

For understanding the concentration of kinetic energy 
along a periodic orbit of characteristics, it is necessary 
to consider in some details the dynamics of these lines. 
Characteristics of PE are, in a meridional plane, straight 
lines making the angle arcsinw with the rotation axis. 
A numerical calculation of their trajectories shows that 
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they generally converge towards a periodic orbit which 
we call, after ]T3), an attractor. The periodic orbit of 
Fig. [| is one example of such an attractor. 




cj = 0.553 
t=-1 



FIG. 1. Kinetic energy in a meridional section of a spherical 
shell of an inertial mode in a viscous fluid. For this numerical 
solution, E = 10 -8 , 570 spherical harmonics and 250 Cheby- 
shev polynomials have been used (the numerical method is 
described in ||). The mode is axisymmetric and symmetric 
with respect to equator. r\ = 0.35 like in the Earth's core, u is 
the frequency of this mode and r its damping rate. Stress-free 
boundary conditions are used. The convergence of character- 
istics towards the attractor is also shown (white lines). 

The Lyapunov exponent (LE) of a trajectory, defined 
by A = liniAr^oo jj Ym=i ln ^ n is the latitude 01 

the n th reflection point), describes how fast characteris- 
tics are attracted or repelled. Its computation as a func- 
tion of frequency shows that attractors (A < 0) are ubiq- 
uitous in frequency space (see Fig. ^). Their existence 
shows that the dynamical system described by the char- 
acteristics is not hamiltonian; the "dissipation" is purely 
geometrical and is due to the fact that, unlike billiards, 
the reflection on boundaries is not specular but conserves 
the angle with the rotation axis. In fact, the dynamics 
of rays is a one-to-one one-dimensional map (from the 
outer boundary to itself), piecewise smooth, but with a 
finite number (twelve) of discontinuities. This kind of 
map has not been studied in the literature of dynami- 
cal systems, perhaps because it does not produce chaos 
because of its invertibility. Iterations of such a map gen- 
erate fixed points which either correspond to attractors 
or to some neutral periodic orbits. Indeed, if r\ = (i.e. 
the sphere is full), all orbits such that to = sm(pn /2q) 
with (p, q) £ IN 2 , are neutral and periodic while those 
such that u> = sin(r7r), r being irrational, are neutral 
ergodic (quasiperiodic). When rj is non-zero only a fi- 
nite number of such neutral periodic orbits subsist; for 
instance, if r\ = 0.35 which is the aspect ratio of the 
Earth's liquid core, q = 1,2,3,4 are the only possibili- 
ties. Interestingly, we face here a situation which is just 



the opposite of the one described by the KAM theorem in 
Hamiltonian systems: when the full sphere is perturbed 
by the introduction of an inner sphere, all ergodic orbits 
are instantaneously destroyed while the longer periodic 
orbits survive the smaller the denominator q is. 




FIG. 2. LE A(cj) of the orbits as a function of u> for 
r\ = 0.35. Inset: blow-up showing the LE of two coexisting 
attractors (full and dashed thick lines). 

Apart from these isolated frequencies which become 
rarer and rarer as r\ increases, generic trajectories are in 
the basin of attraction of attractors. We were able to 
show H that the number of attractors at a given fre- 
quency is finite. The inset of Fig. || shows the typical 
case where an attractor exists in a frequency band [u>x, U2] 
with A(ti>i) = 0, A(cj2) = — 00 and lo2~lo\ ~ 1/N 2 where 
N is the length of the attractor defined as its number of 
reflection points. Near u>x, A ~ ^/lo — u>\ and near u>2, 
A(ui) ~ jj ln N(uj — uo2). The latter implies that long at- 
tractors have small LE in a large fraction of [wi,^] (all 
these results are shown in ffl). 

The existence of attractors for characteristics implies 
that solutions of the inviscid problem (i.e. of PE) are 
singular. This property can be made explicit in the sim- 
plified case of a 2D problem. Indeed, in this case the 
PE may be written d 2 P/du + du- = using characteris- 
tics coordinates; solutions may be constructed explicitly 
from an arbitrary function but, as shown in |HJ, regular 
eigenmodes exist only when neutral periodic orbits exist 
and eigenvalues are infinitely degenerate. When attrac- 
tors are present, the scale of variations of the pressure 
vanishes on the attractors while its amplitude remains 
constant. As velocity depends on the pressure gradient, 
it diverges on the attractor; this divergence is like the 
inverse of the distance to the attractor which makes the 
velocity field not square integrable. This result seems to 
be valid also in 3D j|. 

We therefore understand why solutions of (|l|) look like 
Fig. 1: the inviscid part of the operator focuses energy of 
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the modes thanks to the action of the mapping made by 
characteristics while viscosity opposes to this action via 
diffusion. The resulting picture of Fig. 1 therefore comes 
from a balance between inviscid terms and viscous ones; 
let us make this more quantitative. 
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FIG. 3. Asymptotic behavior of an eigenvalue. The dashed 
line is u — w, as a function of E, while the dotted line is for the 
damping rate r. The solid line represents the 'theoretical' law 
E 1/2 . LUr = 0.403112887 is a root of A(w) = when T) = 0.35. 

For this purpose we first observe that the patterns 
drawn by the kinetic energy of the mode in Fig. [I] is in 
fact a shear layer whose width scale with E a and a ~ 1/4. 
Such a scaling is observed numerically and seems generic 
H[|; it implies that the damping rate of such modes 
scales like E 1 / 2 as clearly shown in Fig. ||. Now we may 
consider a wave packet travelling around an attractor in 
a slightly viscous fluid. The above mentioned balance, 
when applied to both the width and the amplitude of 
the packet, leads to a relation between the LE and the 
Ekman number such as A ~ £ ,1 ~ 3<T with a < 1/3 for an 
eigenmode of the viscous problem. We see that the con- 
straint a < 1/3 is met by actual shear layers. It therefore 
turns out that frequencies of eigenmodes of the viscous 
problem are such that A — * when E — > which means 
that they will gather around the roots of the equation 
A(w) = 0. 

The above result shows the importance of the scal- 
ing verified by shear layers. A boundary layer analysis 
reveals that these shear layers are in fact nested layers 
which consist of an inner a — 1/ 3-layer surrounded by a 
thicker layer. The inner 1/3-layer can be fully explicited. 
Using coordinates along the shear layers (x) and perpen- 
dicular to it (y), we find that the 99-component of the 

velocity verifies ^tjt = — ^t^j with Y = y/E 1 / 3 and 
q = x/y/1 — to 2 which is also the equation verified by 
the stream function in a steady shear layer of a rotat- 
ing fluid fl2]| . Solutions which vanish in Y — ±00 are 
self-similar and of the form u v = q a H a (Y/q 1 / 3 ) with 

H a {t) = / Q °° e- ip *e- p3 p _3a_1 dp. Besides, a = -\ is the 
only admissible value to ensure a coherent evolution of 
the width and amplitudes after reflection on a boundary. 

We are now in position to propose a scenario for the 
asymptotic behaviour of inertial modes when the vis- 
cosity vanishes. Eigenfunctions reduce to nested shear 



layers concentrated along attractors while the associ- 
ated eigenvalues converge toward the frequency u>i such 
that A(u>i) = for the associated attractor. Further- 
more, we can constrain this convergence of eigenfre- 
quencies; indeed, since A ~ y/uj — uji, one finds that 

oj = Ui + aE 2 - 6a H and r = Re(X) = -bE 1 - 2 * when 

E — ► 0; Fig. H shows that this law agrees well with the 
numerical results, in the case shown, with a = 1/4. 

In addition, we noticed earlier that for a finite number 
of to such that ui = sin(p7r/2q) all orbits of characteristics 
are periodic; this implies that in the vicinity of these fre- 
quencies very long attractors with very small average LE 
accumulate as shown by Fig. |j; therefore, these frequen- 
cies will be accumulation points of the asymptotic spec- 
trum. Moreover, around these frequencies eigenmodes 
are weakly damped. On the contrary, modes whose fre- 
quency is in the frequency band of short attractors (like 
the one of Fig. [[]) are more strongly damped. It there- 
fore turns out that the LE curve in Fig. || will strongly 
constrain the distribution of least-damped modes in the 
complex plane at finite viscosities: such modes will avoid 
the large frequency bands of short-period attractors and 
concentrate around frequencies where A(w) = espe- 
cially those with u> = sin(p7r/2q). 

This general evolution of the spectrum is well illus- 
trated in Fig. |[ Here, the least-damped eigenvalues have 
been computed for E = 10~ 8 . We clearly see frequency 
bands of attractors avoided by weakly damped eigenval- 
ues but see the gathering of these eigenvalues around 
sin(7r/4) and, but less conspicuously, around sin(7r/6). 
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FIG. 4. Inverse of the length N of attractors with N < 100 
for r) = 0.35, near the accumulation point 7r/6; each point 
corresponds to an attractor with A = and therefore to a 
point in the asymptotic spectrum. Note the lengthening of 
the attractors as 7r/6 is approached. Here rj = 0.35. 

To complete the picture, we need now mentioning that 
a few regular modes survive among all these singulari- 
ties; such modes are purely toroidal modes or r-modes 
]l3| which are non-axisymmetric. They avoid the con- 
straint of characteristics for their velocity field has no 
radial component; this property makes their characteris- 
tics independent of frequency (they are circles and verti- 
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cal straight lines) and authorizes smooth solution at zero 
viscosity. The associated eigenvalues u> = l/(m + 1), to € 
IN* seem to be the only eigenvalues of the Poincare op- 
erator in a spherical shell. 
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FIG. 5. Distribution of the eigenvalues associated with 
least-damped axisymmetric modes in the complex plane. 
Hatched frequency bands denoted bands occupied by simple 
attractors; the dotted line is for sin(7r/6). The Ekman number 
is 10 -8 and r\ = 0.35. We used a resolution of 700 spherical 
harmonics and 270 radial grid points. 

Ending this letter, it is worth emphasizing the role 
of the geometrical approach allowed by the dynamics of 
characteristics, for describing the asymptotic properties 
of inertial modes; in the domain of very low Ekman num- 
bers (10~ 10 — > 10~ 20 ), typical of astrophysical or geo- 
physical fluids, these modes are out of reach numerically. 

The foregoing presentation shows that inertial modes 
display a very rich dynamical behavior which comes from 
the ill-posedness of the underlying inviscid problem. Here 
we discussed the case of the spherical shell, but our re- 
sults are general and can be extended to any container; 
this is important since natural containers are usually not 
perfect geometrical objects. Hence, fortunately, a curve 
like Fig. || is structurally stable (see our discussion rela- 
tive to the core of the Earth in 

We note that the relevance of attractors has also been 



shown experimentally in stratified fluids |14|. Some con- 
figurations of conducting fluids bathed by a magnetic 
field, obeying the PE, will also display attractors Jl5| . 
These properties are in fact very general and extend to 
mixed-type PDE as illustrated by the case of gravito- 
inertial modes jl6| . We think that similar results should 
hold for systems which are solutions of PDE of hyper- 
bolic or mixed type meeting boundary conditions. As an 
example, our results may have applications in General 



Relativity and the problem of "closed timelike curves" 
(CTC), that is the problem of the existence of physical 
systems which permit travels backward in time. Such 
systems like wormholes have been studied by various au- 
thors p"^ ]; they set many problems among which that of 
causality. Such a problem is also at the origin of the ill- 
posedness of the Poincare problem and we showed that 
it leads to many kinds of singularities. 

We therefore see that inertial oscillations of a fluid in- 
side a container offers a paradigm which may guide our 
intuition for problems in other fields of physics which are 
also ill-posed Cauchy problems. 
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